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We present the impact of two-loop corrections on the mass of the lightest Higgs boson in the 
Minimal R-symmetric Supersymmetric Standard Model (MRSSM). These shift the Higgs boson 
mass up by typically 5 GeV or more. The dominant corrections arise from strong interactions, 
from the gluon and its N = 2 superpartners, the sgluon and Dirac gluino, and these corrections 
further increase with large Dirac gluino mass. The two-loop contributions governed purely by 
Yukawa couplings and the MRSSM A, A parameters are smaller. We also update an earlier analysis 
[1], which showed that the MRSSM can accommodate the measured Higgs and W boson masses. 
Including the two-loop corrections increases the parameter space where the theory prediction agrees 
with the measurement. 


1 Introduction 

The recent discovery at the LHC of a particle consistent with the long sought Higgs boson seemingly 
completes the Standard Model (SM). The mass of the particle is measured with an astonishingly high 
accuracy of run = 125.09 ± 0.24 GeV [2], The precise determination of this mass is of paramount 
importance not only within the context of the Standard Model, but also for finding the path beyond 
it. In fact, a number of experimental observations suggest that the SM cannot be the ultimate theory 
and many theoretical scenarios for the beyond the SM physics (BSM) have been proposed in past 
decades. In some models of BSM, in particular in supersymmetric extensions of the SM, the Higgs 
boson mass can be predicted. However, the current experimental accuracy is far better than theoretical 
predictions for Higgs boson mass in any given model of BSM physics. From the point of view of theory, 
the best accuracy has been achieved in the minimal supersymmetric extension of the SM (MSSM), in 
which the discovery of the Higgs boson and the determination of its mass have given a new impetus 
to the theoretical efforts. The most recent improvements comprise the inclusion of leading three-loop 
corrections [3, 4], resummations of leading logarithms beyond the two-loop level [5, 6], inclusion of the 
external momenta of two-loop self-energies [7, 8], and the evaluation of the 0(af ^contributions in the 
complex MSSM [9, 10]. The MSSM two-loop corrections controlled by Yukawa couplings and a s have 
been known for quite some time for the real MSSM (see the above references for an overview of the 
literature). 

The absence of any direct signal of supersymmetric particle production at the LHC, and the ob¬ 
served Higgs boson mass of ~125 GeV being rather close to the upper value of ~135 GeV achievable 
in the MSSM, are a strong motivation to consider non-minimal SUSY scenarios. In fact, non-minimal 
SUSY models can lift the Higgs boson mass (at the tree level by new F- or D-term contributions 
or at the loop level from additional new states), which makes these models more natural by reduc¬ 
ing fine-tuning. They can also weaken SUSY limits either by predicting compressed spectra, or by 
reducing the expected missing transverse energy, or by reducing production cross-sections. The com¬ 
parison of the measured Higgs boson mass with the theoretically predicted values in any given model 
is therefore highly desirable. Although the theoretical calculations for the SM-like Higgs boson mass 
in such models are less advanced, progress is being made in the development of highly automated tools 
which greatly facilitate the computations in non-minimal SUSY models: SARAH [11—13] automatically 
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generates spectrum generators similar to SPheno [14, 15]; FlexibleSUSY [16] automatically generates 
spectrum generators similar to Softsusy [17]. 

In a recent paper [1] we considered the MRSSM, a highly motivated supersymmetric model with 
continuous R-symmetry [18, 19] distinct from the MSSM. Since R-symmetry forbids soft Majorana 
gaugino masses as well as the higgsino mass term, additional superfields are needed. The MRSSM 
has been constructed in Ref. [20] as a minimal viable model of this type. It contains adjoint chiral 
superfields with R-charge 0 for each gauge sector and two additional Higgs weak iso-doublet superfields 
with R-charge 2. Interestingly, R-symmetry also forbids large contributions to CP- and flavor-violating 
observables [20, 21], so the MRSSM is generically in agreement with flavor data even for an anarchic 
flavor structure in the sfermion sector and for sfermion masses below the TeV scale. Also, Dirac gluinos 
suppress the production cross-section for squarks, making squarks below the TeV scale generically 
compatible with LHC data. Furthermore, models with R-symmetry and/or Dirac gauginos contain 
promising dark matter candidates [22-24], and the collider physics of the extra, non-MSSM-like states 
has been studied [25-33]. 

In Ref. [1] the complete next-to-leading order computation and discussion of the lightest Higgs 
boson and W boson masses has been performed (a similar analysis has been done in Ref. [34]). We 
showed that the model can accommodate measured values of these observables for interesting regions 
of parameter space with stop masses of order 1 TeV. The outcome of the paper was not obvious since 
in the MRSSM (i) the lightest Higgs boson tree-level mass is typically reduced compared to the MSSM 
due to mixing with additional scalars, (ii) the stop mixing is absent and (iii) R-symmetry necessarily 
introduces an SU(2) scalar triplet, which can increase rri\,y already at the tree level. Nevertheless, we 
identified benchmark points BMP1, BMP2 and BMP3 illustrating different viable parameter regions 
for tan /3 = 3, 10, 40 respectively, and also verified that they are not excluded by further experimental 
constraints from Higgs observables, collider and low-energy physics. 

These promising results motivate a more precise computation of the Higgs boson mass in the 
MRSSM and a more precise parameter analysis. Technically, this is facilitated by the Mathematica 
package SARAH, recently updated by providing SPheno routines, which calculate two-loop corrections 
to the CP-even Higgs scalars masses in the effective potential approximation and the gaugeless limit 
[35]. This is the level of precision of the established MSSM predictions except for the refinements 
mentioned above. It is also the level of precision at which the proof [36] applies that the employed 
regularization by dimensional reduction preserves supersymmetry. First applications of the improved 
SARAH version to the calculations of the Higgs boson masses in the R-parity violating MSSM [37] and 
next-to-minimal SSM [38] have been published. 

Since in [1] a judicious choice of the model parameters was needed to meet experimental constraints, 
and an estimate of unknown two-loop contributions was presented, it is of immediate interest to verify 
our findings at higher precision with the new SARAH version. The aim of the current paper is to calculate 
two-loop corrections for the Higgs boson mass in the MRSSM and present an update of the results 
obtained in [1]. 

The paper is organized as follows. After a short recapitulation of the MRSSM in section 2, we 
explain in section 3 our calculation framework and discuss the dependence of two-loop corrections on 
parameters that entered already at the one-loop level. The dependence on parameters that enter only 
at the two-loop level is investigated in section 4. In section 5 we provide an update to the analysis 
presented in [1] using the two-loop corrected masses of Higgses, before concluding in section 6. 

2 The MRSSM 

The MRSSM has been constructed in Ref. [20] as a minimal supersymmetric model with unbroken 
continuous R-symmetry. The superpotential of the model reads as 

IF =/id Rd ■ Hd + Hu Ru ■ H u + A dR-d-T H d + A U R U -T H u 

+ \ d S R d - H d + X U S R u ■ H u -Y d dq ■ H d -Y e el ■ H d +Y u uq ■ H u , (1) 
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where H u d are the MSSM-like Higgs weak iso-doublets, and S, T, R Ut d are the singlet, weak iso-triplet 
and i?.-Higgs weak iso-doublets, respectively. The usual MSSM /j-terrn is forbidden; instead the /J, U ( r 
terrns involving R-Higgs fields are allowed. The A, A-terms are similar to the usual Yukawa terms, 
where the ii-Higgs and S or T play the role of the quark/lepton doublets and singlets. 

The usual soft mass terms of the MSSM scalar fields are allowed just like in the MSSM. In contrast, 
A-terms and soft Majorana gaugino masses are forbidden by R-symmetry. The fermionic components 
of the chiral adjoints, 4>j = O, T, S for each standard model gauge group i = SU(3), SU( 2), 17(1) 
respectively, are paired with standard gauginos g, W, B to build Dirac fermions and the corresponding 
mass terms. The Dirac gaugino masses generated by D-type spurions produce additional terms with 
the auxiliary D-fields in the Lagrangian, 

V n =Mg(B S - V2V b S ) + M^(W a T a - y/2V? v T a ) + M°{g a O a - V2V a g O a ) + h.c., (2) 

which after being eliminated through their equations of motion, lead to the appearance of Dirac masses 
in the scalar sector as well. The soft-breaking scalar mass terms read 


T rEW 
v SB 


m H d (\Hd\ 2 + \Hj\ 2 ) + m 2 Hu (\H°\ 2 + \H+\ 2 ) + [B,(HjH+ - H° d H°J + h.c.] 
+ m 2 Rd (\R°f + \R+\ 2 ) +<X| 2 + m 2 Ru \R~ d \ 2 


+ m 2 s \S\ 2 + m R \T°\ 2 + m R \T | 2 + m R \T + \ 2 + m 2 0 \0\ 2 
+ dL,i m q,i j d L , j + d* R ^m\ i:j d Rij + u* L fn 2 qij u Ljj + u* R ^n 2 Ujij u RJ 
+ ^1^,3 + e%iml i:j e R j + . 


( 3 ) 


The electroweak symmetry breaking (EWSB) is triggered by non-zero vacuum expectation values of 
the R = 0 neutral EW scalars, which are parameterized as 


H d = ^ Vd + ( t >d + ia d) ’ H u = ;j|Kt + &u + i°u) , 

T° = ^( y T + (f>T + icr T ) , s =^(vs + (/>s + icrs)', 

R-Higgs bosons carry R-charge 2 and therefore do not develop vacuum expectation values. We stress 
that in general the mixing of 4>t, <t>s with cj) u and q^ leads to a reduction of the lightest Higgs boson 
mass at the tree-level compared to the MSSM. 


3 Higgs mass dependence on the A, A superpotential parameters 

We now present the MRSSM Higgs boson mass prediction at the two-loop level. We use the same 
renormalization scheme as in Ref. [1], where all SUSY parameters are defined in the DR scheme and 
V S an d vt are determined by minimizing the effective potential at the two-loop order. The 
discussion is divided into two parts. In the present section we begin with the one-loop contributions, 
which are dominated by terms of where a\ collectively denotes squares of the superpotential 

couplings \ u ^ and A u ^- We then discuss the two-loop contributions of 0(a 2 b A ), i.e. ones which depend 
on parameters which already play a role at the one-loop level. In the subsequent section we then discuss 
those two-loop corrections which involve new parameters. 

In the usual MSSM, the one-loop contributions to the Higgs boson mass are dominated by top/stop 
contributions. In the MRSSM, these contributions are also important, but they are simpler since stop 
mixing is forbidden by R-symmetry (corresponding to the MSSM parameter Xt = At — /r/ tan /3 = 0). 
This implies that the top/stop contributions cannot reach values as high as in the MSSM for a given stop 
mass scale. However, as mentioned above, the MRSSM superpotential contains new terms governed 
by A Ut d and A Ut d which have a Yukawa-like structure. References [1, 34] have given a useful analytical 
approximation for these contributions. In the the limit A = X u = —A d, A = A u = Ad, vs ~ v R ~ 0 and 
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Figure 1: Lightest MRSSM Higgs boson mass . and the difference m 2 L ~ miL between masses 
calculated at the two-loop and one-loop level, as a function of A d, X u , A u , respectively. In the upper 
parts of the figure lines from top to bottom correspond to two-loop, one-loop and tree level calculations. 
All other parameters are set to the values of benchmark point BMP1 with tan/3 = 3 (see Tab. 2). 



Figure 2: As in Fig. 1, but for benchmark point BMP3 with tan/3 = 40 (see Tab. 2). 
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This result shows a behavior proportional to A 4 , A 4 and logm 2 oft . This is similar to the top/stop 
contributions as the A’s and 1/ appear in a similar fashion in superpotential. 

We expect therefore that the two-loop result will depend on these model parameters (which already 
entered at the one-loop level) in a manner similar to the pure top quark/squarks two-loop contributions, 
i.e. similar to the MSSM 0(a 2 ) contributions without stop mixing. 

In Figs. 1 and 2 the dependence of the lightest Higgs boson mass calculated at tree-, one- and 
two-loop levels for two benchmarks BMP1 and BMP3 on different model parameters is shown. All 
parameters except the ones shown on the horizontal axes are set to the values of the benchmark points 
defined in Ref. [1] (see Tab. 2). Indeed the A, A behavior of the two-loop corrections is very similar 
to the one of the corresponding one-loop corrections. The numerical impact of the two-loop A, A- 
contributions is rather small, typically less than 1 GeV, except for very large |A U |, |A U | > 1, where they 
can reach several GeV. Particularly, the strong X u dependence for large X u is already manifest for the 
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Figure 3: Feynman rules needed to evaluate diagrams of Fig. 4. In the right diagram, the charge- 
conjugated gluino ^ applies in the case of i = L, gD,0 hi the case of i = R. 


tree-level mass; this is due to the mixing with the singlet state already present in the tree-level mass 
matrix. 

One should remember that very large one-loop contributions are required to bring the predicted 
Higgs boson mass close to the experimental one. In the preferred parameter regions, the A, A are large 
but still moderate enough not to blow up the two-loop contributions. 

Overall, the total two-loop contributions (including the ones to be discussed in the subsequent 
section) are in the range between 4 and 5 GeV, except in the very large A, A regions. This is in 
agreement with the estimate given in Ref. [1], and it confirms the validity of the perturbative expansion 
in spite of the large one-loop corrections. 


4 QCD corrections and the two-loop corrected Higgs boson mass 

At two-loop level the strongly interacting sector and the strong coupling a s appear directly in the 
Higgs boson mass predictions. These two-loop corrections involve not only the gluon but also the 
Dirac gluino and the sgluon, the scalar component of the octet superfield O. They can be expected 
to be sizable, and they depend on the gluino Dirac mass and sgluon soft mass parameters. 1 The 
gluino Dirac mass parameter Mq appears not only directly as the gluino mass but, via Eq. (2), also 
in couplings and mass terms of sgluons, inducing the mass splitting 2 between the real and imaginary 
parts of the sgluon held, O = Os + iOa)- The masses of the scalar sgluons Os and pseudoscalar 
sgluons Oa are related by the tree-level formula tUq s = 4 (Mq ) 2 + m^, w h ere m b A 1S ec l ua ^ to the 
soft-breaking parameter m? 0 [25, 27]. The relevant vertices and Feynman rules are depicted in Fig. 3. 
We assume real Mq , so only the scalar Os acquires the direct coupling to sfermions proportional to 
Mq , via Eq. (2). 

The structure of the strong corrections is thus markedly different from the MSSM case, where only 
the Majorana gluino and the gluon appear. In the following, we study the magnitude and the behavior 
of the corrections as a function of the parameters Mq and m^. 

4.1 Analytic formulas 

As in the previous section, we begin with an analytic approximation for the leading contributions of 
0(octat s ), he. two-loop strong corrections proportional to . This provides us with qualitative insight 
and serves as a check of the code. Generally, in the gaugeless limit the two-loop corrections from 
gluinos and sgluons contribute only to the diagonal part of the {0^, cj) u } submatrix of the scalar Higgs 
boson mass matrix. In the MRSSM the 0(at.a s ) terms contribute only to the (fru&u element. This 

1 These parameters already play a role at lower order, appearing in corrections to Yt (through threshold corrections 
to d s ), though the influence on, for example, DR. top mass is negligible. 

2 In Ref. [1] a simplifying assumption was made that masses of the scalar and pseudoscalar components of (complex) 
sgluon field were equal, since it was unimportant for that analysis. 
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Figure 4: Two-loop diagrams contributing to the Higgs boson mass via Eq. (5) that depend on the 
Dirac mass Mq and the soft sgluon mass mo- We only draw diagrams involving top/stop; similar 
diagrams exist for all quark/squark flavors. 


already constitutes a difference to the MSSM, where the /x-term violates R-symmetry and Peccei-Quinn 
symmetry leading to couplings of stops to < fid- 

Figure 4 shows two-loop diagrams contributing to the Higgs boson mass at 0(ata s ) that explic¬ 
itly depend on mo and/or Mq. These diagrams provide the following contribution to the effective 
potential: 


y( 2 ) _ 

V eff 


D\ 2 


(167T 


2\2 


(MS) 


E 

i=L,R 


fsss(rri~.,m 


> m Os 


) + 


8 gi 


(I 67 P 


? E 

' i=L,R 


fFFs(rnt,m~. 


m 


9d‘ 


(5) 


where the functions fsss and Iffs are dehned in [39]. The effective potential V/ 2 j depends on v u 
through stop masses, which in the gaugeless limit approach 


m lt L 


2 , 1 \/-2 2 
d" 2 *«> 


m inin 


2 1 \z2 2 

m u+2 Y t V u- 


( 6 ) 

(7) 


Equation (5) can be obtained from Ref. [39] by applying translation rules from real fields to complex 
ones. Many such rules can be found in Ref. [35]; an additional rule needed here for the case of a 
Lagrangian C 3 —c<hi|4> 2 | 2 , where 4>i,c G M, 4> 2 G C, is Vsss = \\c\ 2 fsss(m\i "r 2) m^)- 

An important difference to the MSSM is that contributions with fermion mass insertions, corre¬ 
sponding to FFS -type contributions in Ref. [39], are not present in the MRSSM. Such contributions 
vanish due to the lack of L-R mixing between squarks. Hence the gluino mass appears in a simpler way 
than in the MSSM. Likewise, the sgluon only enters via the SSS-type diagram of Fig. 4. An SS -type 
diagram vanishes due to the color structure. 

The corresponding two-loop contribution to the Higgs boson mass matrix element in zero- 

momentum approximation is then given by 3 




(FM _ LAI yf t 

\dv u dv u v u dv u J eff 


( 8 ) 


For large tan/3, corrections of order 0(at,a s ) cannot be neglected any more. But since they contribute 
only to (frd&d matrix element, their impact on mass of the lightest Higgs, which stems mainly from the 
4> u (t>u element, is small. Results of Eq. (5) where compared with the results of two-loop routines from 
the SARAH-generated SPheno module. 
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Figure 5: Two-loop contributions to the SM-like Higgs boson mass depending on the gluino mass in 
the MRSSM for BMP1 (left panel) and BMP3 (middle panel) and two different values of the soft 
sgluon mass parameter mo = 2 TeV (thick solid blue line) and 10 TeV (thick dashed green line) 
with all contributions, respectively, and without the sgluon contributions (thin solid red line). For 
comparison also the MSSM contributions for no (thin dashed light blue line) and maximal (purple 
dotted line) stop mixing are plotted. The chosen MSSM parameters are given in Tab. 1. For BMP3 
and mo = 2 TeV the right panel shows the result, when successively switching off dominating and 
sub-dominating contributions. 

_| tan/3 Mi M 2 /r m A m^ M;(3|3) m lu,d m je A r,b X t 

left panel 3 600 500 400 700 1000 2 2000 2 1000 2 0 0/2000 

right panel 40 250 500 400 700 1000 2 2000 2 1000 2 0 0/2000 


Table 1: Definition of the fixed parameters for the MSSM points in Fig. 5. All parameters in GeV or 
GeV 2 , where appropriate. The stop mixing parameter Xt is given both for the cases of no and large 
stop mixing. 

4.2 Numerical analysis 

We now turn to the numerical analysis of the complete two-loop corrections to the SM-like Higgs boson 
mass, using the full evaluation within the framework of SARAH and SPheno. The first two panels of Fig. 5 
focus on the gluino and sgluon mass dependence, which arises mainly from the 0(ata s ) corrections; 
they show the two-loop corrections as a function of the gluino mass parameter for two different values 
of the soft sgluon mass, mo = 2 and 10 TeV for two benchmarks BMP1 and BMP3; other parameters 
are fixed at benchmark values. For comparison, the two-loop result without the sgluon contribution is 
shown as well (i.e. without the first diagram of Fig. 4). We also plot the MSSM prediction with strong 
stop mixing and without any sfermion mixing the at tree-level. 

The first two panels show that the dependence in the MRSSM without sgluon contributions is 
very similar to the one in the MSSM without stop mixing. The corresponding thin solid red and 
thin dashed light blue curves in Fig. 5 show a characteristic drop for large gluino masses. This is 
understandable as in the MSSM without sfermion mixing the gluino contribution is precisely the same 
as in the MRSSM and given by the two corresponding diagrams in Fig. 4. The Dirac or Majorana 
nature of the gluino does not matter since the Dirac partner, the octet superfield O has no direct 
couplings to quark superfields. A few TeV gluino masses slightly increase the Higgs boson mass, but 
for larger values of Mq the /ffs function becomes negative and drives the correction downwards. 

In the full MRSSM calculations, including the sgluon diagrams strongly changes the behavior. 

3 As pointed out in [35], in SARAH and SPheno the two-loop tadpole contributions are included directly in vacuum 
minimization condition and not in Eq. 8. 
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Surprisingly, the full MRSSM two-loop contributions resemble the MSSM contributions with large 
stop mixing. In both cases, large gluino masses strongly enhance the Higgs boson mass, however for 
different reasons. In the MSSM the increase can be traced back to the additional FFS -type diagram 
which is directly proportional to Mq and which vanishes in the limit of no stop-mixing. In the MRSSM, 
on the other hand, the sgluon diagram grows with Mq both due to the sgluon-stop-stop coupling, which 
scales like Mq , and to an increase in the scalar (but not pseudoscalar) sgluon mass. Due to the sgluon 
contributions the total two-loop contributions to the Higgs boson mass in the MRSSM are larger than 
the ones in the MSSM. They are further increased by heavy sgluons. 

The third panel of Figure 5 compares the numerical impact of individual contributions by succes¬ 
sively switching off contributions. It allows to read off the contributions from sgluon, gluino and gluon, 
of 0(af,atab), and the remaining two-loop contributions (particularly the A, A contributions). The 
gluon diagrams alone contribute approximately +4 GeV. The negative gluino and the positive sgluon 
corrections together amount to an additional upward shift of the Higgs boson mass, which can reach 
several GeV for large Dirac gluino masses. The remaining contributions are far smaller and amount to 
around —1 GeV for the 0(a^,atCXb) contributions and +0.5 GeV for the remaining contributions. 


5 Update of benchmarks 
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Figure 6: Contour plots showing the behavior of rriH 1 given by the color map and mw given by the red 
contour lines. The plots are ordered horizontally by benchmark points (top row for BMP1, bottom row 
for BMP3), while vertically different combinations of model parameters are varied. The white stars 
mark the original benchmark points from Ref. [1], whereas the black ones show the adapted points 
after taking into account the two-loop corrections. 

In this section we present an update of the analysis of Ref. [1], using the more precise evaluation 
of the Higgs boson mass. Ref. [1] studied the mass predictions of the W and lightest Higgs bosons in 
the MRSSM and showed that agreement with experimental data is possible, in spite of tree-level shifts 
from violations of custodial symmetry and from mixing with other Higgs states, respectively. 

Table 2 shows benchmark parameter points defined in that reference. They exemplify parameter 
regions in which mw and rri //, agree with experiment. They are characterized by large |A| ~ 1, rather 
light Dirac higgsinos and gauginos, and they have tan (5 = 3,10,40, respectively. 

For all three benchmark points the two-loop correction to rri is around +5 GeV. As discussed in 
the previous sections, the largest part of this is due to the 0{ata s ) corrections. The MRSSM-specific 


8 





































BMP1 

BMP2 

BMP3 

tan (3 

3 

10 

40 

B, 

500 2 

300 2 

200 2 

^u 

1.0,-0.8 

1.1,-1.1 

0.15,-0.15 

A-di A u 

-1.0,-1.2 

-1.0,-1.0 

-1.0,-1.15 

Mg 

600 

1000 

250 

2 

rrir> 

flu 

2000 2 
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Iddi 
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K 
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Mg 
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m^p. m|, m 2 0 

3000 2 ,2000 2 ,1000 2 

m Q;l,2> m Q ;3 

2500 2 ,1000 2 

™>D-, 1,2> m D ;3 

2500 2 ,1000 2 

m U -,3 

2500 2 ,1000 2 

m'i: m| 

1000 2 

2 

Trim 

fid 

700 2 


v s 

4.96 

0.67 

-0.30 

Vt 

-0.34 

-0.20 

-0.34 

m H d 

673 2 

743 2 

1160 2 

2 

TTlrr 

flu 

—535 2 

—542 2 

—541 2 

m H i 

130.3 GeV 

130.3 GeV 

129.8 GeV 

mw 

80.400 GeV 

80.384 GeV 

80.393 GeV 

HiggsBounds’s obsratio 

0.67 

0.68 

0.67 

HiggsSignals’s p-value 

0.03 

0.03 

0.03 


Table 2: Benchmark points of Ref. [1], Dimensionful parameters are given in GeV or GeV 2 , as appro¬ 
priate. The first two parts define input parameters. The third part shows parameters derived from 
electroweak symmetry breaking after solving the tadpole equations at two loops. The last part gives 
the theory predictions for the Higgs boson mass at the two-loop level and further quantities relevant 
for comparison with experiment. 


corrections of are small since the values of A u , though large, are still not as large as needed to 

make these corrections dominate, see Figs. 1, 2 for two out of three benchmarks. The magnitude of 
the total two-loop correction is consistent with the theory error estimate given in Ref. [1], 

The upward shift of rriH 1 implies that it is easier to obtain agreement with the measured value, 
i.e. smaller values of |A U | are sufficient. In Tab. 3 we provide new, slightly modified benchmark 
points, whose definitions differ only in the values of A u . The two-loop Higgs boson mass prediction 
agrees well with experiment, and the good agreement of mw with experiment is unchanged. Likewise, 
both the old and the new set of benchmark points pass checks against HiggsBounds [40-42] and 
HiggsSignals [43, 44], 

In Fig. 6 we give an update to some of the subfigures from Figs. 4 and 5 of Ref. [1], These show the 
predictions of mw and mn 1 as contour lines in several two-dimensional parameter spaces. The Higgs 
boson mass is evaluated at the two-loop level. As discussed before, with the exception of the regions of 
very large A, there is a general positive contribution to the lightest Higgs boson mass between 4 and 
5 GeV. Accordingly, the contour lines, in particular the central green region in which the Higgs boson 
mass agrees with experiment, shift to slightly lower values of A. Also, the overlap region, where Higgs 
and W boson masses agree with experiment, is enlarged. 
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BMPT 

BMP2’ 

BMP3’ 

A-U 

-1.11 

-0.85 

-1.03 

V S 

5.2 

1.01 

-0.22 

Vt 

-0.25 

-0.02 

-0.21 

m H d 

674 2 

764 2 

1160 2 

mir 

£1 u 

—502 2 

—512 2 

—516 2 

m Hl 

125.3 GeV 

125.5 GeV 

125.4 GeV 

m\v 

80.397 GeV 

80.381 GeV 

80.386 GeV 

HiggsBounds’s obsratio 

0.61 

0.65 

0.87 

HiggsSignals’s p-value 

0.72 

0.66 

0.72 


Table 3: Adapted benchmark points; other parameters are as given in Tab. 2 


6 Conclusions 


In this work we have presented the impact of two-loop corrections on the mass of the lightest Higgs 
boson in the MRSSM. The calculation has been performed using the framework of SARAH in the 
approximation of the vanishing electroweak gauge couplings and external momenta of the Higgs self 
energies. The code has been cross-checked with an analytic calculation of the most important new 
corrections. We have separately analyzed the impact of contributions involving the A, A-couplings, 
which already appear in the one-loop corrections, and of the strong corrections involving gluon, Dirac 
gluino, and sgluon exchange. 

In the previous work [1] and the present paper we have found that the lightest Higgs boson mass 
in the MRSSM differs from the one in the usual MSSM in several respects. At tree-level the additional 
mixing with additional scalar states reduces the MRSSM Higgs mass below the MSSM value. At the 
one-loop level, the top/stop contributions cannot be as large as in the MSSM, because stop mixing is 
forbidden by R-symmetry. However, the new contributions from the superpotential A, A-terms have a 
similar structure as the top/stop contributions. If the A, A-couplings are similar in magnitude to the 
top Yukawa coupling, the lightest Higgs boson mass can easily be in the ballpark of the experimentally 
allowed range. 

The two-loop corrections governed by these A, A-couplings, however, amount to only 1 GeV or less 
in parameter regions in which the Higgs boson mass agrees with experiment. The most important 
two-loop contributions are the strong corrections of 0(ata s ). As we have shown the Dirac gluino and 
gluon contributions alone are very similar to the MSSM strong contributions for vanishing stop mixing. 
The inclusion of the sgluons changes the picture. The sgluon contributions are positive and rise with 
the Dirac gluino mass, such that the total 0(ata s ) corrections of the MRSSM are larger than the ones 
of the MSSM, independently of the magnitude of stop mixing. 

Overall, the MRSSM two-loop corrections to the lightest Higgs boson mass are typically positive. 
E.g. for the benchmark parameter points proposed in Ref. [1], the two-loop corrections to the Higgs 
boson mass amount to approximately +5 GeV, within the error estimate of that reference. Since 
perturbation theory shows a converging behavior and since the A, A-corrections are subdominant (for 
|A|, |A| less than around 1.2), we estimate the remaining theory uncertainty to be not larger than the 
one of the MSSM. 

The positive two-loop corrections make it easier to achieve agreement between the theory prediction 
for the lightest Higgs boson mass and the measured value. We have provided an update of the analysis 
of Ref. [1], showing parameter regions of simultaneous agreement of the Higgs and W boson mass 
predictions with experiment. Compared to Ref. [1], the allowed parameter regions are slightly larger 
and located at smaller values of the A, A-couplings. 
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